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Abstract 

We derive the equations for approximate particle number projection based on mean 
field wave functions with finite range density dependent forces. As an application 
ground bands of even-A superdeformed nuclei in the A = 150 and A = 190 re- 
gions are calculated with the Gogny force. We discuss nuclear properties such as 
quadrupole moments, moments of inertia and quasiparticle spectra, among others, 
as a function of the angular momentum. We obtain a good overall description. 
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1 Introduction 



Understanding pairing correlations has been and still remains a challenge for 
nuclear physicists. Several decades ago it was thought that pairing correlations 
will only affect a few nuclei and the properties of a few observables at low 
energy and small angular momentum; mainly the level density of odd nuclei 
as compared with the even ones, the energy of the two quasiparticle states and 
the moments of inertia of rotational nuclei. Such simplistic view has changed 
over the years in such a way that nowadays we expect pairing correlations 
playing a major role in many nuclear phenomena. To mention a few recent 
examples let us remember the important role of pairing correlations in halo 
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nuclei, in N = Z systems and in the understanding of the second moment of 
inertia in superdeformed nuclei. 

In spite of all that the theoretical understanding of the pairing force is rather 
limited. The recognition of this fact has motivated in the last years a prolifer- 
ation of approximate methods (Lipkin-Nogami, Kamlah at second order and 
exact particle number projection) and new pairing forces (density dependent, 
surface and volume pairing). 

Within the mean field approaches the most successful calculations, aiming at 
describing low energy nuclear structure phenomena have been performed with 
density dependent forces. The most popular, because of its simplicity, is the 
Skyrme force. Most of the calculations have been performed with the SkM* 
parameterization [1,2]. The main disadvantage of the Skyrme interaction is 
that realistic pairing matrix elements cannot be extracted from the force, and 
therefore the pairing interaction (normally a monopole pairing term) has to 
be added by hand. Recently, other alternative forms of the pairing interac- 
tion, like density dependent pairing have been explored in connection with 
the Skyrme interaction [3-5]. Theories beyond mean field with simple pairing 
forces have been also developed [6]. Pairing properties, on the other hand, are 
very important in the description of high spin properties. Since the pioneer 
work of Mottelson and Valatin [8] many papers and experiments have been 
devoted to this topic. This is why the investigation of high spin states with 
a realistic force that is able to treat the short range correlations properly is 
most desirable. The Gogny force [9] is well known for providing, with the same 
set of parameters, a good description of many nuclear properties at low spins 
over the nuclide chart as, for example, ground state energies, odd-even mass 
differences, electron scattering data [10], fission barriers [11], etc. The main in- 
gredients of the force are the phenomenological density dependent term which 
was introduced to simulate the effect of a G-matrix interaction and the finite 
range of the force which provides the pairing matrix elements from the force 
itself. In recent calculations the Gogny force has been applied to high spin 
problems at normal deformations [12] and to superdeformed shapes [13,14]. 

In these calculations the cranked HFB approximation (CHFB) was used. The 
angular momentum and number of particles of the system were determined 
by the corresponding constraints on the wave function. It turned out [12], 
however, that in spite of the finite range of the force, the condition of holding 
the particle number to be right on the average was not enough and a superfluid 
normalfluid sharp phase transition was found, contrary to what one would 
expect for a finite system. To remedy this problem, in [15] we have formulated 
the Lipkin-Nogami method for density dependent forces, the so-called CHFB 
Lipkin-Nogami approximation(CHFBLN). When this formalism was applied 
the sharp phase transition was washed out, the results of the calculations being 
in much better agreement with the experiment. 
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The purpose of this paper is to report on calculations in the A = 190 and 
A = 150 regions within the CHFB and CHFBLN approximations and the 
Gogny force. In section 2.1 we first derive the Lipkin-Nogami approximation 
from the Kamlah expansion to second order. In sect. 2.2 we extend this for- 
malism to density dependent forces obtaining an extra contribution due to the 
rearrangement term. In sect. 3.1 we discuss results for the A = 190 region, 
namely for the 19 °- 192 - 194 Hg and 194 Pb, and in sect. 3.2 the A = 150 region, 
taking as examples the 152 Dy and 150 Gd nuclei. We conclude the paper with 
conclusions and some appendixes. 



2 Projected Mean Field Theories 

Mean field theories (MFT) are based on the assumption of a product wave 
function; they constitute the basic (and in most cases the only one feasible) 
approximation to a many body system. The success of MFT is based on the 
concept of spontaneous symmetry breaking. By this mechanism the Hilbert 
space is considerably enlarged and many correlations can be taken into ac- 
count. Obviously, the corresponding wave functions have to be considered as 
"intrinsic" wave functions. The most sophisticated mean field theory is the 
HFB approximation, in which the intrinsic wave function | $) is defined as the 
non-trivial vacuum to the quasi-particle annihilation operators 

^ = E^ + ^c+, (i) 



where |cj,c^~j is the set of single particle creation and annihilation operators 
in terms of which the quasi-particle operators are expanded, U and V are the 
so-called Bogoliubov wave functions. Product wave functions, |$), of quasi- 
particles ct M , in general, are neither eigenfunction of the angular momentum 
nor of the particle number operator. For heavy, well deformed nuclei, the usual 
approach of demanding that, at least, the expectation values of the angular 
momentum and the particle number operators take the right value, is a good 
approximation for the angular momentum but not for the particle number. 
This has to do with the number of nucleons participating in the breaking of 
the associated symmetry [17]. 

The ideal treatment of pairing correlations in nuclei is particle number pro- 
jection (PNP) before the variation [17]. At high spin this theory is rather 
complicated to implement and up to now it has only been applied to separa- 
ble forces [18]. On the other hand, the semi-classical recipe of solving the mean 
field equations with a constraint on the corresponding symmetry operator can 
be derived as the first order approach to a full quantum-mechanical expansion 
(the Kamlah expansion [19]) of the projected quantities. The second order in 
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this expansion takes into account the particle number fluctuations and might 
cure some of the deficiencies of the first order approximation. However, full 
selfconsistent calculations up to second order are rather cumbersome and just 
simple model calculations have been carried out up to now [20]. Most sec- 
ond order calculations have been done using the Lipkin-Nogami (LN) recipe, 
which was originally proposed in Refs. [21-23]. Recently, the LN method has 
been applied to study superdeformed nuclei at high spins adding a monopole 
pairing interaction to the Woods-Saxon plus Strutinsky approximation [24], 
and to the Skyrme force [25]. 

The usual derivation of the LN method is not valid for density dependent 
forces because the dependence in density is not taken explicitly into account. 
In refs. [15,16] the LN approach was generalized to properly deal with this 
situation. The derivation is based in the Kamlah expansion for the particle 
number projected energy. In the new formulation additional terms arise in 
the equations determining the parameters of the theory. To introduce the 
approximation and the corresponding notation we make in the following a 
short derivation of it. 



2.1 First and second order approximation to exact particle number projec- 
tion : Derivation for non-density dependent forces 

The LN method can be derived using different arguments. As the usual ways 
of deriving the LN method are not well suited for a generalization to the case 
of density dependent forces, in this section we shall derive the equations for 
non-density dependent forces using the Kamlah expansion. This derivation 
can be easily extended to the density dependent case, see sec. 2.2. 

Let |$) be a product wave function of the Hartree-Fock-Bogoliubov type, i.e. 
a particle number symmetry violating wave function. We can generate an 
eigenstate of the particle number, \^n), by the projection technique [17] 





o 



The particle number projected energy is given by 



E: 



'proj 



(V N \H\y N ) _ / o 2 "d0e-^ JV ($|Fe^|$) _ Jo" d<Pe- i4,N h(<p) 
Wn\^n) ~ It #e-^($|e^|$) ~ Jt d<Pe-^ N n((j)) 



, (3) 
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where we have introduced the Hamiltonian and norm overlaps 

h((f)) = ($|#e^|$), n(0) = ($|e^|$). (4) 

In the case of large particle numbers and strong deformations in the gauge 
space associated to N, one expects the h((p) and n((f>) overlaps to be peaked 
at = and to be very small elsewhere in such a way that the quotient 
h((f)) / n(4>) behaves smoothly. In this case one can make an expansion of h(<f>) 
in terms of n(0) in the following way [19] 

M 

h(<j>) = ]T h m Ar m n(<l>), (5) 

m=0 

where we have introduced the Kamlah operator 

tf=±JL-(q>\N\Q), (6) 

which is a representation of the particle number operator in the space of the 
parameter 0. M represents the order of the expansion, for M — * oo eq. (5) 
is exact. For large symmetry breaking systems and large particle number one 
expects convergence with a few terms. The expansion coefficients h m are de- 
termined by applying the operators 1, N, J\f M on eq. (5) and taking the 
limit — > 0, one obtains 

M 

(<f>\H(AN) n \$) = h m (^\(AN) m+n \^) (7) 

m=0 

with AN = N — (<&| jV|<&) and n — 0, . . . , M. The previous equality provides 
an M + 1 equations system, ho, hi, ... , hu being the unknowns. From now on 
we will use the shorthand notation (A) = Substitution of eq. (5) in 

eq. (3) provides the projected energy to order M 

M 

Kroj = E h m (N - (N)r (8) 

m=0 

This expression tell us that to determine the wave function |$), in an ap- 
proximate (to order M) projection before the variation, we have to minimize 
eq. (8). To first order, i.e. M — 1, one obtains 

= h + hi(N-{N)), (9) 
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with 



(10) 



In this case to minimize the approximate projected energy, we have to solve 
the HFB equations with the constraint (N) = N, i.e. the old BCS recipe. This 
approximation is expected to be good for very strong violating wave functions. 
In atomic nuclei where only a few pairs participate in the pairing correlations, 
we do not expect the first order approximation to be a good one. 

For M = 2 one obtains the following equations system for the coefficients 
ho, hi and h 2 

(H) = h +0 + h 2 ((AN) 2 ) 

< (HAN) = + h((AN) 2 ) + h 2 ((AN) 3 ) (11) 

(H(AN) 2 ) = h ((AN) 2 ) +h 1 ((ANf) + h 2 ((AN)*) 

the solution of this equations system provides 

h =(H)-h 2 ((AN) 2 ) 
_ (HAN)-h 2 ((ANf) 
((AN) 2 ) 

_ ((H- (H))(AN) 2 ) - (HAN)((ANf)/((AN) 2 ) 
2 ((AiV) 4 ) - ((AiV) 2 ) 2 - ((AN) 3 ) 2 /((AN) 2 ) 

The projected energy to second order is given by 

E% = (H) - h 2 ((AN) 2 ) + h x (N-(N)) + h 2 (N 

In a full variation after projection method one should vary E^} j and compute 
cumbersome expressions like f|. In the Lipkin-No garni prescription, however, 
the coefficient h 2 is held constant during the variation. As a result the varia- 
tional equation is much simpler, one gets 

4<#> " ^i-(N) + (N — (N)) 5 ^ - 2h 2 (N - (N))4-(N) - h 2 ^-((AN) 2 ) 



(12) 
(13) 

(14) 



-<A0) 2 (15) 



It is obvious that the solution to this equation is given by 



±(H - h 2 (AN) 2 ) - \±(N) = 0, (17) 



6 



with A determined by the constraint 

(N) = N , 



(18) 



provided the condition A = hi is accomplished. This can be easily realized 
noticing that eq. (17) must hold for any variation |<5$), in particular, we can 
choose \5<&) = AiV|$). Substitution of this specific variation in eq. (17) pro- 
vides 

(HAN}-\{(AN) 2 }-h 2 ((AN) 3 } = (19) 
comparison with eq. (13) shows that A = hi. 

We have formulated the approximation for one class of nucleons (i.e. protons 
or neutrons), in practical cases one has to do it simultaneously for protons and 
neutrons but it is easy to show that for M = 2 there is no coupling between 
protons and neutrons. The energy to be minimized is 

E$ oj = (H) - hl({AZf) - ^((AN) 2 ), (20) 

with the constraints (Z) = Z and (N) = N. In these expressions Z stands for 
protons and N for neutrons, h\ is given by eq. (14), where the number oper- 
ator in this equations should be substituted by the proton number operator, 
h 2 is defined correspondingly. In the same manner if there are additional con- 
straints, for example the angular momentum, one just has to substitute H by 
H' — H — uJ x in all equations above with u chosen in such a way that the 

2 

constraint (J x ) = [1(1 + 1) — (J z )] 1 ^ 2 is accomplished. 

2.2 First and second order approximation to exact particle number projec- 
tion: Derivation for density dependent forces 

Now we would like to generalize the formulae of the former section to density 
dependent forces, like the Gogny force [9-11] 

vn(p) = E e^ 1 -^ 2 ^ (Wi + B t P a - H t P T - MiP P T ) 

i=l,2 

+ iW LS (V 12 5(fi - f 2 ) A Vi 2 )(ai + <t 2 ) 

+ t 3 (l + P a x )6(r 1 -r 2 )p 1 / 3 ( i ^p), (21) 

where the dependence on the wave function |$), through the density p(r) = 
($|ct(f)c(f)|$), adds some peculiarities. 
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Density dependent forces do not exhibit, in general, the same properties as 
the bare two-nucleon force. Let |$) be a ^-symmetry breaking wave function 
and V\ 2 (p) the corresponding interaction. It is obvious that 

Sv 12 (p)S' 1 = v 12 (S P S- 1 ). (22) 

That means, the two-body interaction in a S-rotated system is the same as 
the interaction calculated with the rotated density. Eq. (22) is valid for any 
symmetry S, in the particular case of symmetries with no dependence on the 
spatial coordinate r, for example the particle number, it also holds 

SvMS- 1 = v^SpS- 1 ) = v 12 (p), (23) 

i.e. the Vi 2 (p) interaction and the generators of the symmetry S do commute. 
In particular, for the particle number, in which we are interested here, V\ 2 (p) 
satisfies the same properties as the bare force. 

In the context of plain mean field theories there is no ambiguity left concerning 
the density dependence of the hamiltonian. In theories beyond mean field, for 
example the Generator Coordinate Method or projected theories, in evaluating 
the expectation value of the energy, appear matrix elements like (&\H\&), or 
in the case of PNP, the overlap ($\He l< t' N \$). Here a complication arises from 
the density dependent term of the hamiltonian because it is not obvious which 
density dependence should be taken in the potential (21). A Skyrme force with 
a linear dependence in the density can be cast as a three body force. For this 
particular case it can be shown that the density to be used in the hamiltonian 
in the evaluation of the matrix elements (&\H\&) must be a non-diagonal 
density Pnd(j) given by 

pNo{r) = = ($|$') » (24) 

where we have introduced the density operator p(r) = c + (r)c(r). For the par- 
ticular case of PNP, this density is given by pnd{t) = ($\p(r)e i4>N \®) / (<&\e i<t>N \<&) 
Up to now, based on this result, the same prescription has been adopted for 
Skyrme forces with a non linear dependence in p [26] as well as with the Gogny 
force [15]. This prescription is, however, not well founded. In the context of 
the LN approximation there have been two approaches to this problem. The 
first one [15], considers the energy as a functional of the densities (normal 
and abnormal) while the second one [16], is based on a many-body density 
dependence of the potential. These two methods are discussed in the following 
subsections. 
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2.2.1 Functionals of the density 



In this section we shall derive the approximate particle number projection 
equations for potentials with a density dependence like eq. (21), using the 
prescription of appendix B to calculate the non-diagonal matrix elements of 
the hamiltonian. There the use of a density like (24) in the evaluation of 
non-diagonal matrix elements of any density dependent Hamiltonian was mo- 
tivated. 

The density dependence of the hamiltonian does not modify the projected 
energy of eq. (8) but only the coefficients h m of the expansion. The non- 
diagonal density causes a dependence on <fi of the hamiltonian and eq. (5) 
looks like 



M 



($|#(0)e^|$> = Y. h m N m (<5>\e^ N \<5>), 



(25) 



m=0 



application of the operators 1,/v, ...,Af M on eq. (25), in the limit — > 0, now 
provides 



E 

i=o 



n 



1 d 1l ~ l ~ 



M 



) = £ h m ({AN) m+n ). (26) 

0=0/ m=0 



For M — 1 the energy is given by eq. (9) and the coefficients ho and hi by 



IdH 

i d<p 



(H) = h 

+ HAN\ = h 1 {(AN) 2 }, 
=o / 



(27) 



the derivatives with respect to are calculated very easily, since 



ldp 

i d(j) 



1 d ($|p(f)e^|$) 



0=o * d( t> ($|e^|$) 



= (p(r)AN) 



(28) 



we obtain 



IdH 

i d(p 



IdH dp 

=0 i dp dtp 



(29) 
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Substitution in eq. (27), provides the coefficient hi 



($\HAN\<S>) 
($|(AA>) 2 |$) 



dH(p) ($|pAiV|$) 



dp(r) ($|(AiV) 2 |$) 



(30) 



This expression for hi is exactly the same obtained in the selfconsistent HFB 
approximation, see eq. (61). It is interesting to notice that if we had used the 
density ($|p(r|$), instead of eq. (24), in the hamiltonian, the coefficient hi 
would be given by the first term of eq. (30), i.e. the rearrangement term would 
be missing, contrary to the HFB result, see appendix A. 

For the M = 2 case, we obtain 



1 d 2 H 



i 2 d<j) 2 



lldH 
ldH 



(H) = h + h 2 {(AN) 2 ) 
+ HAN ) = hi({AN) 2 ) + h 2 ((ANf 



+ 2AN . 

=o 1 d $ 



(31) 



+ H{AN) 2 ) = h ((AN) 2 ) + hdAN) 3 ) + h 2 ((AN) 4 



0=0 



The second derivative with respect to is given by 

1 d 2 p 



i 2 d(f) 2 



(Ap(ANY) 



0=0 



(32) 



and 



1 d 2 H 



i 2 d<j) 2 



(33) 



From now on we proceed as in the non-density dependent case, i.e. we have 
to solve eq. ( 17 ) with the constraint (18) but with the coefficient h 2 given 
by the solution of equations system (31). It can be easily checked that the 
condition X — hi also holds in this case. 



2.2.2 Many-body density dependence of the Hamiltonian 

In the former section we have interpreted the Gogny potential as a functional 
of the density and we have derived, under this assumption, the approximate 
particle number projected equations. 

Another approach [16] to the same problem is to look at the expectation value 
(^jvl-f^l^jv) of eq. (3) independently of the fact that \^n) is a many body wave 
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function. In this case one could think of using the density (^n\p\^n) /(^n\^n) 
in the hamiltonian. This density, which we shall refer to as projected density 
is given by 

_ {<!>\p n pwp»\$) _ m P (riP N m 



and should be used in the Gogny force (21)0 in the calculation of the pro- 
jected energy. We shall call Hj(= H(p pro j)) the hamiltonian obtained with a 
dependence in the projected density. We would like to remark that the wave 
function |$) in eq. (34) is a symmetry breaking one, i.e. it is not an eigenstate 
of N. This guarantees the use of generalized Hamiltonians which are function- 
als of wave functions in the very large Hilbert space generated by symmetry 
breaking wave functions. Concerning the projection, we are now able to calcu- 
late the projected energy of eq. (3) and we can proceed along the same lines 
as in the derivation for non-density dependent forces, i.e. to use the Kam- 
lah expansion up to order M. All expressions are the same, one just has to 
substitute H by Hi in all equations. 

Calculation of the projected density may be quite involved. A way out is to 
approximate the projected density to be used in the hamiltonian. This can be 
done in the following way. Substitution of (2) in p pro j{f) provides 



PproAT> / o 2 "d0 e -^($|e^|$) " Jo" d(j)e-^ N n((j)) ' 1 ' 

where we have introduced the overlap p{r,4>) = ($|p(r)e JI?uv |$). As with the 
Hamilton overlap, we may use the Kamlah expansion to approximate the 
projected density, i.e. 

p(f,<P)= E R m {^ m <(t>). (36) 

m=0,M 



Substitution of eq. (36) in eq. (35) provides the projected nuclear density to 
order M 

M 

&>E^(^-W) m (37) 

m=0 

2 This prescription, which, as we shall see works for the particle number, may not 
work for symmetries depending on r. This has to do with the fact that the particle 
number symmetry is a special one because the spatial density is invariant against 
rotations in the gauge space associated with N, see eq. (23). 
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The expansion coefficients R m are determined applying the operators 1, A/, A/ M 
on eq. (36). In the limit 0^0, one obtains 

M 

($|/3(f)(AiV)"|<f) = £ iUr)($|(AiVr + "|$>, (38) 

m=0 



with n — 0,...,M. To first order 



= ^o(r) + i2i(f)(JV - (A>», (39) 



with R = ($|p(r)|$) = p(f), the HFB density, and 



V ; ((AN) 2 ) V ; 



Assuming that R\(r)(N — (N)) is a small quantity as compared to p(r) we 
can expand the density which appears in the Gogny force 



[p$ oj (r)} 1/3 «P(0 1/S 



p(r) 



1/3 



«p(r') 1/3 + ip(f)- 2 / 3 i2i(f)(JV- <A>», (41) 

substitution of this expression in the Gogny hamiltonian with the projected 
density dependence, what we have called Hj, provides 

Hj = H( Pproj ) » (p) + — i?!(f)(iV - (A>», (42) 



with -ff(p) the hamiltonian depending on the HFB density. The projected 
energy is given, up to (N — (N)) order, by 



E% = (Hj) + ^^(iV - (AO) = (H) + hf(N - (N)) (43) 



with hf 1 = ^ivp) + (^f-^i(O)- Minimization of this energy is equivalent to 
the minimization of (H) — X(N) subject to the constraint (N) = N and A 
given by h\^. Now A (or h\^) includes the rearrangement term present in 
the HFB theory with density dependent forces. It is worth to point out that 
in this approach we obtain the same results as in the previous section , see 
eq. (30), and in the selfconsistent HFB approximation, see eq. (61). 
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To calculate the energy to second order one should also expand the density 
overlap to the same order. We obtain 

= P{r) + Ri(r)(N - (N)) + R 2 (r)(N - (N)) 2 - it> 2 (f)((AiV) 2 )(44) 



with 



= (p^AN) - R 2 (r) ((ANf) 
U ' ((AiV) 2 ) 1 ' 

R (f) = (W) - (p(r1))(AiV) 2 ) - (p(r)AN)((ANf)/((ANr) 

((AiV) 4 ) - ((AiV) 2 ) 2 - ((AiV) 3 ) 2 /((AiV) 2 ) 1 ' 



As in eq. (41), we can make an expansion of the dependence of the hamiltonian 
on the projected density. After substitution of this expansion in the projected 
hamiltonian, we obtain 

Hj » H(p) + [ Rl (f)(N - (N)) + R 2 (f)(N - (N)) 2 - R 2 (f)((AN) 2 }} (47) 



The projected energy for M = 2 is given by 

E%L = (Hi) hi((AN) 2 } + h[ (N — (N)) + h{ (N - (iV)) 2 (48) 



with 



(W)-/4((A^) 
((AiV) 2 ) 

h i = m - (H^jAN) 2 } - (^AiV)((AiV) 3 )/((AiV) 2 ) 
2 ((AiV) 4 ) - ((AiV) 2 ) 2 - ((AiV) 3 ) 2 /((AiV) 2 ) ' 1 ' 

The structure of eq. (48), is exactly the same as eq. (15) and so are the 
coefficients h[, h\ and h±, h 2 - From the discussion of eq. (15), it is obvious that 
the solution to this equation is given by 

^<iJ 7 -/4(AiV) 2 )-A^(iV) = 0, (51) 



with A determined by the constraint 

(N) = N, (52) 
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provided the condition A = h[ is accomplished. Again, it can be checked that 
this is the case. Taking into account the particle number constraint we can 
write 

Hi^H{p)-^^R 2 {r){{^Nf)) 



Substitution of Hj in h 2 provides 



hI = h yfej^g) ~ (|fei? 2 (f)))(AiV) 2 ) - (^iZ 2 (f)AJNr)((AJV)3)/((AJV)") 
2 2 ((AiV) 4 ) - ((AiV) 2 ) 2 - ((AiV) 3 ) 2 /((AiV) 2 ) 

= h 2 + <ta 2 (53) 

with h 2 given by eq. (14). We have checked in some test cases that 8h 2 is 
much smaller than h 2 and therefore in the following 5h 2 is neglected. Assuming 
(N) = N, we can finally write the projected energy, eq. (48), to be minimized 



f)FT 

E% « (H) - K-Q^Mrl) + h 2 ]((ANY) 

= (H)-h e 2 ff ((AN) 2 ). (54) 
Where we have introduced the effective parameter h e 2 ^ 



h e ff = ((H-(H))my) - (HAN)((ANf)/((AN)*) 
2 ((AiV) 4 ) - ((AiV) 2 ) 2 - ((AiV) 3 )7((AA0 2 ) 



with 



K^H + ^i—Mr))^ (56) 



The quantities fij(r) are those appearing in the second quantization form of 
the density operator p = Ylij fij(^) c t c j- The last term in this equation is 
a consequence of the density dependence of the hamiltonian and does not 
appear in the standard derivation of the Lipkin-Nogami approximation. It 
resembles the usual rearrangement term appearing in the HFB theory with 
density dependent forces. This term might be very important for forces with 
density dependent pairing. 
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3 Results for ground bands of even-even nuclei. 

In the evaluation of the LN energy the most complicated expression is (^AH(AN) 2 
in appendix C we derive the general expression for this term. The detailed ex- 
pressions for the total LN energy and its gradient for the Gogny force can 
be found in ref. [27]. The equations have been solved using the Conjugate 
Gradient Method [28]. 

We have applied the LN formalism with the Gogny force to study the superde- 
formed ground bands of even A nuclei. We shall refer to this calculations as 
cranked-HFB-LN (CHFBLN). In most results we shall also present the ones 
obtained with the plain cranked HFB theory (CHFB) [12]. We are using the 
standard D1S parameterization set [11,29]. In some figures the theoretical re- 
sults and the experimental data are displayed as a function of the angular 
momentum, in spite of the fact that for some nuclei the spins have not been 
firmly assigned. The experimental spin values that appear in these figures are 
those quoted in the respective works as the most probable ones. 

To solve the CHFB and CHFBLN equations we have expanded the quasipar- 
ticles in a triaxial harmonic oscillator (HO) basis. The HO configuration space 
was determined by the condition 

huj x n x + hoOyUy + hu z n z < hu N (57) 
where n x , n y and n z are the HO quantum numbers and the frequencies hu x , hu y 

1 2 

and hu z are determined by uj x = uo y = Uoqs, cu z = u q~z. The parameter q is 
strongly connected to the ratio between the nuclear size along the z- and the 
perpendicular direction. A value of q = 1.5 is a good election for superdeformed 
nuclei. For N we have taken a value of N = 12.5, which provides a basis big 
enough as to warrant the convergence of the results. 

As it is usually done, to save CPU time, the following terms of the interaction 
have not been taken into account in the calculations : The Fock term of the 
Coulomb interaction and the contributions to the pairing field of the spin- 
orbit, Coulomb and the two-body center of mass correction terms. 

It is important to notice that in the calculations of the coefficient hi of the 
Kamlah expansion, the full hamiltonian is used in the expressions. Concerning 
the two prescriptions to deal with the density dependence of the hamiltonian, 
in ref. [16] it was shown that for the A = 190 region both prescriptions provide 
very similar results. Therefore, for this region we show results only with the 
prescription of the many body density dependence of the hamiltonian. For the 
A = 150 region we shall show results for both prescriptions. 
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3.1 The A = 190 region. 



In this section we shall present results for the ground band of the nuclei 
190-192-194 anc j 194P]-, j n ]_ we ^qt^ the particle-particle correlation 

energy, defined by E pp = |7V(A/c), in the CHFB and in the CHFBLN ap- 
proaches, versus the angular momentum, for the nuclei 19 °- 192 - 194 Hg. 




8 16 24 32 40 48 8 16 24 32 40 48 8 16 24 32 40 48 
1(H) 1(h) 1(h) 

Fig. 1. Particle-particle correlation energies for 19 °-i92-i94jjg^ -p u jj triangles : CHFB 
pairing energy for protons (upwards) and neutrons (downwards). Empty triangles : 
CHFBLN pairing energy for protons (upwards) and neutrons (downwards). 

In the CHFB approach we do not find pairing correlations for the proton sys- 
tem for the nuclei considered. For the neutron system we find sizable pairing 
correlations at low and medium spins which diminish at large angular mo- 
mentum due to a strong Coriolis antipairing effect, at very high spins they get 
very small. In the CHFBLN approach, we find larger correlation energy for 
both, the proton and neutron systems, we also observe the Coriolis antipairing 
effect as the angular momentum increases but the systems remain correlated 
even at the largest spins. We also realize that the LN term has a larger effect 
on the proton system than on the neutron one. 

The energy E pp is a measure of the pairing correlations, but there are other 
ways to look for pair correlations. In the scientific literature, and in the scope 
of HFB theories, the pairing energy Ep is defined as 

Ep = E — Enp (58) 

with E the energy calculated in the HFB (or LN) approximation and E HF the 
energy calculated in the Hartree-Fock approach. This quantity is different from 
what we have called particle-particle correlation energy E pp . The quantities 
Ep and E pp coincide only in the case of a BCS calculation performed on top of 
a HF solution, i.e. without selfconsistency in the Hartree-Fock and the pairing 
fields. 
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In fig. 2, as an example, we display the quantity Ep for 190 Hg, in the HFB 
and the LN approach for protons and neutrons. In the HFB solution we have 
found a superfluid solution in the neutron channel but not in the proton 
channel. The HF solution provides the proton and neutron densities, pz and 
Pn, that minimize the energy in the absence of pairing correlations. We expect, 
therefore, the proton contribution to the total energy to be deeper in the HF 
than in the HFB approach, the neutron contribution, on the contrary, must 
be deeper in the HFB than on the HF approach. The total energy must be, 
obviously, deeper in the HFB than in the HF approach. We observe in fig. 2 
that these expectations are fulfilled. We also find a slow convergence to zero 
in the energy Ep, as a function of the angular momentum, indicating the 
quenching of the pairing correlations. In the LN approach we find the neutron 
system to behave very similarly to the HFB one, the proton one behaves 
differently as expected from the behavior of E pp . It is interesting to point out 
that E P , a quantity with more physical content than E pp , is much smaller 
than E pp as one would expect from the experimental data. It is important to 
notice that whereas the correlation energy gained by the HFB procedure with 
respect to HF is only of 0.5 MeV at / = 12h, the correlation energy gained 
by the LN procedure with respect to HFB at the same spin is of 2.0 MeV. 
This result indicates the relevance of the restoration of the particle number 
symmetry in this case. 

In the upper part of fig. 3 we show the angular momentum dependence of the 
quadrupole moment, Qo p , for the nuclei 19 °- 192 - 194 Hg. With increasing neutron 
number we observe a global decrease in deformation and a reduction in the 
/-values in which the quadrupole moment starts to decrease. The CHFBLN 
and the CHFB results differ mostly at low spins where the pairing correlations 




12 16 20 24 28 32 36 40 44 48 

1(H) 

Fig. 2. Pairing energies, eq. (58), for 190 Hg. Full triangles : CHFB pairing energy for 
protons (upwards) and neutrons (downwards). Empty triangles : CHFBLN pairing 
energy for protons (upwards) and neutrons (downwards). 
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12 24 36 48 12 24 36 48 12 24 36 48 
I (h) I (h) I (fi) 

Fig. 3. Results for !90-i92-i94jjg Tjpp er p ane l : The charge quadrupole moment 
versus the angular moment. Full (empty) circles represent the CHFB (CHFBLN) 
approach. Lower panel : Gamma-ray energy AEj = E(I) — E(I — 2) as a function of 
the angular momentum. Full (open) circles stand for the CHFB (CHFBLN) results 
and full squares for the experimental data, for 190 Hg [31], for 190 Hg [33] and for 
194 Hg [34] . 



are larger. This causes the crossing of the results of both calculations for each 
nucleus. As a function of the angular momentum we observe first an increase of 
the quadrupole moment due to a reduction of the pairing correlations and then 
a decrease of Qq p . This anti-stretching effect is caused by the Coriolis force. The 
CHFB and the CHFBLN predictions are slightly higher than the experimental 
values 17.7li. 2 eb for 190 Hg [30], 17.7 ±0.8 eb for 192 Hg [32] and 17.7 ±0.4 eb 
for 194 Hg [32] . These quadrupole moments correspond to (3 deformations along 
the yrast band of 0.560 - 0.580, 0.555 - 0.565, 0.540 - 0.552 for i 9 o-i92-i94 Hg 
respectively, the gamma deformation remains zero for the three nuclei at all 
spin values. 

In the lower part of fig. 3 we display the transition energies versus the angular 
momentum. The agreement with the experimental data is good in the CHF- 
BLN approach, in the CHFB one obtains too low values. The reason for this 
behavior is again the smallness of the pairing correlation in the CHFB ap- 
proach as it can be seen in the upper part of fig. 4 where we display the static 
moments of inertia versus the square of the angular frequency Hcu. Here we 
observe that the CHFB moments of inertia are too large as compared with the 
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Fig. 4. Upper panel : the static moment of inertia, 2^\l) = (2I—l)/AEj versus the 
square of the angular frequency, Twj(I) = AEj/2. Lower panel : Second moments of 
inertia, l (2 \l) = 4/(A£j- A£ 7 _ 2 ), for i90-i92-i94 Hg _ The meaning f the sym bols 
is the following : Squares for the experimental data [31], [33], [34], full (empty) circles 
for the CHFB (CHFBLN) approach. 



experiment. The CHFBLN moments of inertia are in much better agreement 
with the experiment, specially for 192 - 194 Hg. 

In the lower part of fig. 4, the second moment of inertia X^> as a function 
of the angular frequency is represented. Here, as in the previous picture, the 
CHFB results are in poor agreement with the experiment. At small spin values 
the moment of inertia is too large and the neutron alignment sets in too early. 
The CHFBLN results, however, are in very good agreement with the exper- 
imental results, the correlations induced by the Lipkin-Nogami prescription 
produce the desired effect, namely, to diminish the moment of inertia at small 
spin values and to delay the alignment to higher spin values. For 190 Hg the 
experimental results show a sharp upbending at / = A2h ( hu ~ 0.4 MeV ), 
in our calculations we do not get a sharp upbending but a modest one, we 
see however a clear discontinuity in the 1^ at I = AAh. For the last points 
of our calculations there are not experimental results. One expects, however, 
that after the sharp backbending the experimental results will bend down. 
The CHFBLN for 192 Hg are in very good agreement with the experimental 
ones and for 194 Hg we get slightly more alignment at high spins than in the 
experiment. 
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0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

Tico (MeV) Tico (MeV) 

Fig. 5. Quasiparticles energies for the nuclei 190 Hg (upper panel) 192 Hg (middle 
panel) and 194 Hg (lower panel), in the left hand side for positive parity protons and 
in the right hand side for negative parity neutrons. The continuous (dashed) lines 
correspond to positive (negative) signature. 

To investigate the orbitals responsible for the alignments we display in fig. 5 the 
quasiparticle energies for the three nuclei. We find that the orbitals responsible 
for the alignment at Hou ~ 0.4 MeV are the N = 7 ones, mainly the i/[761]| for 
190 Hg and 192 Hg and the z/[752]| for 194 Hg, though it is hard to distinguish both 
orbitals because they are strongly mixed. These results are in good agreement 
with the ones obtained in [3]. The proton alignment takes place at higher 
frequencies (huj ks 0.45). 

As a further example of this region we display in fig. 6 the results for the 
nucleus 194 Pb in the CHFBLN approximation. The particle-particle correla- 
tion energy E pp , fig. 6a, shows a similar behavior to the Hg isotopes, i.e., a 
modest Coriolis antipairing effect at small spins which becomes larger at high 
angular momentum. The quadrupole moment behaves similarly to the previ- 
ous studied nuclei : first, it increases as a consequence of the weakening of 
the pairing correlations and at high angular momentum starts to decrease due 
to the anti-stretching effect. The numerical value of Q 0p is also a little larger 
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Fig. 6. The nucleus 194 Pb in the CHFBLN approximation, (a) The particle-particle 
correlation energy for protons (triangles) and neutrons (inverted triangles) (b) The 
charge quadrupole moment, (c) the transition energy, the experimental data are 
from [36] and (d) the dynamical moment of inertia 

than the experimental value of 18.8 ±1.1 eb, [35]. The transition energies are 
depicted in fig. (6c), here the agreement with the experimental data is excel- 
lent. Finally, in fig. (6d) the dynamical moment of inertia is shown, we find 
a very good agreement with the experimental data at low and high angular 
momentum. The analysis of the quasiparticles energies as a function of the 
angular momentum, indicates that in this case the orbital responsible for the 
alignment at high spin is the z/[752]|. 



3.2 The A = 150 region. 



We shall now turn to the A = 150 region. This region is characterized by 
smaller pairing correlations, the restoration of particle number symmetry 
should not be, therefore, as relevant as in the A = 190 region. As examples of 
this region we discuss the nuclei 152 Dy and 150 Gd. 

As mentioned at the beginning of this section, we would like to compare, in the 
A = 150 region, the results obtained with the two approaches to the problem 
of the density dependence of the Hamiltonian discussed in the theory section. 

In fig. 7 we represent the main results for the nucleus 152 Dy in the HFB and 
in the LN approach. The LN results with the prescription of the many-body 
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Fig. 7. The nucleus 152 Dy in the LN and the HFB approximations, (a) The par- 
ticle-particle correlation energy in CHFBLN (CHFBLNFD) for protons, open tri- 
angles (circles) and neutrons, empty inverted triangles (squares). Same quantity in 
CHFB : protons, triangles; neutrons, inverted triangles, (b) The charge quadrupole 
moment, CHFBLN (empty circles), CHFBLNFD (squares) and CHFB (full circles), 
(c) the transition energy. Experiment (full squares) [41], CHFBLN (empty circles), 
CHFBLNFD (empty squares), CHFB (full circles) and (d) the dynamical moment 
of inertia. Same symbols as in (c). 

density dependence we shall call CHFBLN, as before, and those based in the 
functional of the density prescription, CHFBLNFD. 

The particle-particle correlation energy, E pp , is represented in fig. 7a as a 
function of the angular momentum. In the CHFB approximation, the neutron 
system is not correlated while the proton one is slightly correlated for spin 
values smaller than 36h. In the LN approximation, however both systems are 
correlated. The results in the CHFBLNFD and CHFBLN are very similar for 
protons and neutrons. As expected we obtain smaller correlation energies in 
this region than in the A = 190 one. In fig. 8 we plot the energy E P of eq. (58) 
in the CHFBLN approach. For neutrons is very small for all spin values and 
for protons around 1.8 MeV at spin 40h and it decreases to 1.5 MeV for spin 
60h, i.e. the energy gain allowing superfluid wave functions is about 2 MeV. 

In fig. 7b we display the charge quadrupole moment, Q 0p , versus the angular 
momentum. In the CHFB approximation and at small angular momentum it 
is almost constant, only around spins values of 30h, Qo p starts to decrease by 
the Coriolis anti-stretching effect. In the CHFBLN approaches the behavior at 
low spins is quite different (it increases due to the Coriolis antipairing effect) 
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Fig. 8. Pairing energies, Eq. (58), for 152 Dy in the CHFBLN approximation for 
protons (dots) and neutrons (squares). 

and only around spin values of 50h starts to decrease. The experimental value 
is 17.5 ± 0.2 eb, [39,40], in good agreement with the CHFBLN results, the 
HFB values being somewhat larger. The nuclei are axially symmetric, i.e. the 
gamma deformation is zero. In fig. 7c the gamma ray energies are plotted 
versus the angular momentum. The HFB predictions are the lowest ones due 
to the large moment of inertia. The CHFBLN results are slightly below the 
experimental data and almost parallel to them. Both CHFBLN approaches 
to the gamma ray energies practically coincide. Finally in fig. 7d the second 
moments of inertia are displayed versus the angular frequency. The CHFB 
results are larger than the experimental ones but they behave in a similar 
way. The inclusion of additional pairing correlations through the LN approach 
produces an enhancement of the second moment of inertia as compared with 
the CHFB results. This anomalous behavior has been also found in [7] . The two 
LN approaches are very similar at high spin values, only at small spins one finds 
differences. The low values of X*- 2 -* at small angular frequencies are probably 
due to an excess of pair correlations at low /, as one can also see in the behavior 
of the quadrupole moment, fig. 7b, being a second derivative magnifies 
small effects. The first moment of inertia, fig. 7c, however, has very reasonable 
values. Concerning the quality of the LN approaches, we may conclude, as in 
the A = 190 region that both approaches practically coincide. 

In fig. (9) the quasiparticle energies for the nucleus 152 Dy are represented 
versus de angular frequency. From these levels one can predict the blocking 
structure of the excited bands in this nucleus or in the nearby odd-even nuclei. 
One should be aware, however, that there is strong mixing within the Nilsson 
levels and that the character of particle or hole near the Fermi surface, in the 
presence of pairing correlations, is not a precise statement. Most of the levels 
shown in the figure coincide with the assignment proposed [41,42]. 
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Fig. 9. Quasiparticles energies for the nucleus 152 Dy, in the upper part for protons 
and in the lower part for neutrons. The continuous (dashed) lines correspond to 
positive (negative) signature. 

The next example in the A = 150 region is the 150 Gd, the results are shown in 
fig. 10. Since we have seen in the previous example that both density prescrip- 
tions provide similar results, for this nucleus we only show the results with 
the many body density prescription. In fig. 10a, the CHFBLN pairing energies 
versus the angular momentum are represented, they show an smooth decrease 
with growing angular momentum. In fig. 10b, the charge quadrupole moments 
as a function of / are presented. They display the typical Coriolis effects: first 
an increase caused by the Coriolis antipairing effect and later one the Corio- 
lis antistretching effect. It is in good agreement with the experimental value, 
eb, [37]. This nucleus is also axially symmetric. The transition energies 
represented in fig. 10c are well reproduced by the theoretical predictions, indi- 
cating, nevertheless, a slightly larger moment of inertia than the experiment. 
The second moment of inertia is shown in fig. lOd, the agreement is good ex- 
cept at the lowest values of the angular momentum, where the experimental 
points present a sharper level crossing than in the theory. 



4 Conclusions 



In conclusion, we have discussed the Lipkin-Nogami equations for density de- 
pendent forces and two different prescriptions for the density dependence of 
the Hamiltonian to be used in the evaluation of non-diagonal matrix elements. 

As an application we have analyzed the high spin properties of the nuclei 
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Fig. 10. The nucleus 150 Gd in the CHFBLN approximation, (a) The particle-particle 
correlation energy for protons (triangles) and neutrons (inverted triangles) (b) The 
charge quadrupole moment, (c) the transition energy, the experimental data [38] 
are represented by squares and (d) the dynamical moment of inertia 



i90-i92-i94jjg an( j i94pk Q f A = 190 region and the nuclei 152 Dy and 
150 Gd of the A = 150 region in the Cranked Hartree-Fock-Bogoliubov and the 
Cranked Lipkin-Nogami approaches. 

Global properties as well as sensitive properties as the second moment of in- 
ertia are well reproduced, specially within the CHFBLN approximation. In 
general we obtain excellent qualitative agreement with the experiment and 
in some cases a quantitative one. The approximate restoration of the particle 
number symmetry is the essential ingredient for the quantitative agreement. 
We believe that the Gogny force with the approximate particle number pro- 
jected mean field theory provides a good description of nuclear properties 
under a large variety of situations. 

This work was supported in part by DGICyT, Spain under Project PB97- 
0023. One of us (A.V.) would like to thank the Spanish Ministerio de Asuntos 
Exteriores for financial support through an ICI grant. 
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5 Appendix A :Lagrange multiplier with density dependent forces 



In the HFB theory the wave function |$) that minimizes the energy with the 
constraint of providing the right particle number on the average is given by 

5{<Z>\H - A7V|$) = (59) 



for density dependent forces one obtains 

dH(p) 



(5$\H\<S>) + ($ 



dp(r) 



mm 



$) - A(5$|iV|$) = 



(60) 



this equation must hold for any arbitrary variation |<5$), in particular for 
|<5$) = AiV|$), substitution in eq. (60) provides 



A 



($\HAN\$) 
($|(AiV) 2 |$) 



+ $ 



dH(p) ($|pAJV|$) 



dp(r) ($|(AiV) 2 |$) 



(61) 



6 Appendix B : Functionals of the density 

We want to define a prescription to evaluate overlaps of the hamiltonian 

for density dependent forces where both |<3>) and |<E>') are product wave func- 
tions of the HFB type. The philosophy behind the prescription is that a density 
dependent force does not define a hamiltonian but rather the energy, which is 
a functional of the density and the pairing tensor 

E [Pi K \ = T klk2 p k2kl + - v k 1 k 2 h A k A (2pfe 4 fc 2 Pfc 3 fei + K k 3 k 4 ^t lk2 ) (62) 

fcl&!2 klk2k$k4 



where the two-body interaction could depend upon the spatial density p(r) 
to some power a, v[p a (f)}, as is the case for the Gogny force. In the above 
expression the spatial density is defined in the usual way 

p( r ) = ~J^y = 2- /*ib( r >b*i' 
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where the density operator is given by 

A 

pif) = Y / s(f- fi) = /fcife 2 (^)4 1 cfc 2 - 

i=l kik 2 



The formal evaluation of the hamiltonian overlap for non-density dependent 
forces is straightforward. Using the extended Wick theorem, one gets 

^ k\k 2 kzk/t {^Pk 4 k 2 Pk i k 1 + Kk 3 k 4 K kik 2 ) ' 



$\H\& 

\ I / k\k 2 k\k 2 kzk 4 



where 



Therefore, for non-density dependent forces both, the energy, eq. (62) and the 
hamiltonian overlap, eq. (63), are given by the same functional of the density 
matrix and the pairing tensor but evaluated with the corresponding densities 
and pairing tensors. Inspired by these results we use the following prescription 
for density dependent forces : in the evaluation of the hamiltonian overlap a 
density dependent interaction depending upon p, eq. (64), must be used. From 
the point of view of defining a hamiltonian the previous prescription amounts 
to use the density 



for the density dependent part of the hamiltonian. This prescription coincides 
with the one obtained for a Skyrme force with a linear dependence in p. 
The above density is, in general, not a real quantity rendering the density 
dependent "hamiltonian" a non hermitian operator. This is not a real problem 
as the hamiltonian overlap can be a complex quantity, we only have to be sure 
that the projected energy is a real quantity as it will be shown below. 

Another important check concerning the consistency of the prescription is the 
following : since P N and the density dependent hamiltonian do commute and 

N ~ (MM ~ ' { } 
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According to our prescription, to evaluate non-diagonal matrix elements of 
the first expression of in eq. (65), in the hamiltonian we have to consider 
the density 



($|e-^p(f)e^ ,iV |$> 
($| e iW-4>)N |$) 



(66) 



whereas for the non-diagonal matrix elements of the second expression of 
in eq. (65) one should use 



The final result, should be independent of whether we use the first or the 
second expression. This is clearly true, since eq. (66) can be written as 



due to the commutativity of e l ^ N with /5(f). The operator e l ^ N also commutes 
with po^'-^r) showing that 



h(</>, 0') = ($| e~^H[p^,(r)}e^* |<t>> = ($| H[p^, ^{r)^' ~^ |$> 
= MO,0'-0) 



which is the property needed to obtain the second expression of Eq (65). Obvi- 
ously, we shall use the second expression for the projected energy and eq. (67), 
for the density entering in the density dependent part of the interaction, in 
the evaluation of h(<f>). 

Next, in order to ensure the reality of the projected energy, we have to show 
that the following property holds 




(67) 



Po^'-^(r) 



(Q\p{r)jW-® N \Q) 



h*{4>) = h{-4>). 



(68) 



The complex conjugate of p^if) is given by 




($| e-^ N p{r) |$) _ {<$>\ p{r)e-^ N |$) 
($| e -^|$) " ($| e -^|$) 
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using this property in the definition of h*(4>), Eq (68) is demonstrated, 
h*{4>) = ($| e-^H{]?Ar)) |$) = ($| H{p^{f))e-^ |d>> = h(-< 



7 Appendix C: Evaluation of (AH(AN) 



In order to evaluate (AH (AN) 2 ) we will use the following relation 



AH(AN) 2 ) = 









e i(pN 


5>> 



(69) 



(f>=0 



In this expression the density appearing in the density dependent part of the 
hamiltonian is, obviously (<3>|c + (f)c(r)|$). The matrix element of the right 
hand side is given by eq. (63) with |$') = e % * N \&). We may further write 



($|e^ |$) 



Tr 



-Tr 



A(0)k(0) 



where we have introduced the notation 

ks k4 kg 



Performing the second derivative with respect to <fi and taking the limit = 
we obtain 



d 2 ($|ife^|$) 



= -Tr 



Tp 2 N + ^2NP + Fp 2 N + ^NPn) 



+ -Tr[A 2N K* + A7t 2JV + 2A n k n ], 



where we have introduced the matrices 



(70) 



(pN) klk2 ~ ~ ^Pfcxfc 2 (0) 

(«Jv)fc 1 fc 2 = 7 



</>=o 



0=0 



( /c ^)fcifc a = 7 ^"fcitaW 



(cj>| c+c fcl A7V|$>, 
($| c fc2 c fel A7V|cJ>>, 
($| C +c+A7V|<J>>, 
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and 



(<t>| A (c+c fel ) (AiV) 2 |$), 

($|A(c fc2 c fel )(A7V) 2 |$), 

($|A( C +c+) (AiV) 2 |$), 

rAT,r 2 Ar, Atv and A 2 tv are given by the usual expressions but the traces have 
to be taken with Pn,P2N, k n and k 2 n respectively. 
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